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Abstract: We study the spectral behavior of higher order elliptic operators upon 
domain perturbation. We prove general spectral stability results for Dirichlet, Neu¬ 
mann and intermediate boundary conditions. Moreover, we consider the case of the 
bi-harmonic operator with those intermediate boundary conditions which appears in 
study of hinged plates. In this case, we analyze the spectral behavior when the bound¬ 
ary of the domain is subject to a periodic oscillatory perturbation. We will show that 
there is a critical oscillatory behavior and the limit problem depends on whether we 
are above, below or just sitting on this critical value. In particular, in the critical case 
we identify the strange term appearing in the limiting boundary conditions by using 
the unfolding method from homogenization theory. 
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1 Introduction 

In this paper, we consider the general problem of the spectral behavior of an elliptic 
partial differential operator (i.e., the behavior of its eigenvalues, eigenfunctions as well 
as of the solutions to the corresponding Poisson problem) when the underlying domain 
is perturbed. In with A > 2, we will consider a family of domains {D£}o<e<eQ 
which approach a limiting domain D as e —)• 0, in certain sense to be specified and we 
will also consider higher order selfadjoint operators (order 2m with m > 1) with not 
necessarily constant coefficients and with certain boundary conditions. The operators 
will have compact resolvent and therefore the spectrum will consist only of eigenvalues 
of finite multiplicity. 

Importantly, the associated energy spaces, generically denoted by V(De), will satisfy 
the condition C V{Qe) C and will depend on the domain and 

the boundary conditions considered. We will consider different types of boundary 
conditions according to the choice of the spaces V(De). If V(De) = they 

will be called Dirichlet boundary conditions, if P(De) = IT™'’^(De) they will be called 
Neumann boundary conditions, and in case V{Q^) = n IPQ’^(De) for certain 

1 < k < m — 1, they will be called “intermediate boundary condition”. We refer to 
[7] and references therein for a pioneer discussion on the stability properties under the 
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three types of boundary conditions, including an analysis of the so-called Babuska- 
Sapondzhyan paradox. We also refer to [29] for a further discussion on the paradox 
and |30] for a general reference in this type of problems. We mention that sharp 
stability estimates for the eigenvalues of higher order operators subject to Dirichlet 
and Neumann boundary conditions have been recently proved in m where uniform 
classes of domain perturbations have been considered (see also m for related results); 
moreover, in mm further restrictions on the classes of open sets allow obtaining also 
analyticity results. 

Our goal is twofold. On one hand, we will provide a rather general condition 
describing the way the domains converge to the limiting one, which will guarantee 
the spectral convergence of the operator in Og to the appropriate limiting operator 
in 0. The condition, which we will denote by (C), see Section [3] below, is expressed 
intrinsically and it is posed independently of the boundary conditions that we consider. 
Needless to say that for a particular family of perturbed domains to check that the 
condition is satisfied will depend heavily on the boundary condition imposed. This 
condition generalizes previous ones formulated for Dirichlet and Neumannn boundary 
conditions. 

On the other hand, we will focus on the case of higher order operators with “inter¬ 
mediate boundary” conditions, paying special attention to the case of the biharmonic 
operator. We will obtain almost sharp conditions on the way the boundaries can be 
perturbed to guarantee the spectral convergence with preservation of the same inter¬ 
mediate boundary conditions for general higher order operators. Afterwards we will 
analyze in detail the case of the biharmonic when the boundary of the domain presents 
an oscillatory behavior. We will see that there is a critical oscillatory behavior such 
that, when the oscillations are below this threshold we have spectral stability, while 
for oscillations above this value we approach a problem with Dirichlet boundary con¬ 
dition. For exactly the threshold value, there appears an extra term in the boundary 
condition for the limiting problem, which maybe interpreted as a “strange curvature”. 
The existence of this critical value is well known in other situations. See for instance 
the seminal paper [20] and also, |29] . In other context see HaElEH]. 

We describe now the contents of the paper. In Section [2] we set up the operators, 
fix the notation and include a subsection where we describe the basic elements of 
the technique called “compact convergence of operators” which will be used in this 
paper. In Section [3] we state condition (C), see Definition 13.II and show that this 
condition implies the “compact convergence of operators” and therefore, the spectral 
stability of the operators. In Section [4] we consider the case of Dirichlet boundary 
conditions while in Section [5] the case of Neumann boundary conditions. The case of 
intermediate boundary conditions is studied in Section [6l We prove Lemma 16.21 and 
Corollary 16.181 which provides conditions guaranteeing the spectral stability for the 
intermediate boundary conditions. These sections cover the first goal of the paper. 

The second goal is achieved in sections [7] and [8] . Notice that in Section [7] we analyze 
the case of a biharmonic operator with intermediate boundary conditions (in this case 
F(D) = IF^’^(D) n IFq’^(D)) where the domain is perturbed in an oscillatory way. As 
a matter of fact if the boundary of the unperturbed domain is given locally around 
certain point by the function xn = ^(aJi, • • • ^xn-i) for x = (xi,... ,xn-i) £ W for 
some nice (iV —l)-dimensional domain W, with 5 = 0 (that is, the boundary is flat) and 
the boundary of the perturbed domain is given as xn = ge{x) where 5 e(x) = e"6(x/e) 
for some smooth and periodic function 6, then a = 3/2 is a critical value. If a > 3/2, 
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the oscillations are not too strong and the limit problem has also the same intermediate 
boundary conditions. If a < 3/2, the oscillations are too wild and the limit problem 
has a Dirichlet boundary condition in W. The critical case a = 3/2 is treated in detail 
in Section [8j We need to treat this case as a homogenization problem and will use the 
unfolding operator method to show that the limit boundary condition in W contains 
an extra term. The results of this paper were announced in [5]. 


2 Preliminaries, notation and some examples 

We introduce in this section the general setting of the paper, the basic notation to 
follow the contents and some relevant examples. Also, we include the definition of 
“compact convergence”. 


2.1 Higher order elliptic operators 

We fix some notation and we recall basic facts from standard spectral theory for elliptic 
operators. We refer to Davies [22] for details and proofs. 

Let N,m G N and Q be an open set in M^. We denote by the Sobolev 

space of real-valued functions in L^(D), which have distributional derivatives of order 
m in L^(D), endowed with the norm 

— II^IIl2(o) + ^ (2.1) 

\cx\=m 

We denote by the closure in 1T™'’^(D) of the space of the C'°°-functions with 

compact support in D. 

Let rh be the number of the multi-indices a = (ai,..., oat) G Nq with length 
|a| = ai ajsf equal to m. Here No = N U {0}. For all Nq such that 

l®l = 1/^1 = '^1 ^a/3 be bounded measurable real-valued functions defined on 

satisfying = AjSa and the condition 


Aa/3(x)^aCl3 > 0 

\a\ = \P\=m 


for all X € ^ = { ia )\ a\=m £ For all open sets D in we define 

Qn{u,v) = ^ / AapD°‘uD^v dx + / uv dx, 

, , ^ Jn Jn 


\a\ = \P\=m' 


( 2 . 2 ) 


(2.3) 


for all u, u G VF™'’^(D) and 

Qn{u) = Qn{u,u). 

Observe that by condition (12.21) Qq is in fact a scalar product in ^"^’^(D). 

Let F(D) be a linear subspace of IT™’^(D) containing IT((”’^(D). We recall that if 

1 /2 

F(D) endowed with the norm (•) is complete then there exists a uniquely deter- 

1 /2 

mined non-negative selfadjoint operator such that Dom Hyp^ = F(D) and 

1/9 1/9 

Quiuiv) < Hy^^^u > 1,2 p'j , V'u,uGH(D). 
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In particular, a function u belongs to the domain of if and only if tt € 1^(fl) and 

there exists / G such that 


Qn{u,v) =< f,v VuGl/(fI). (2.4) 

Clearly, Hy^Q'^u = /. Equation (12.4p is the weak formulation of the classical problem 

Lu = /, in n (2.5) 

where L is the classical operator defined as 

Ln:=(-1)™ +u (2.6) 

\a\=\l3\=m 


and the unknown u is subject to suitable homogeneous boundary conditions depending 


on the choice of E(n) (see the examples below). 

We recall that if the embedding E(n) C is compact then the operator 

has compact resolvent. In this case the spectrum is discrete and consists of a sequence 
of eigenvalues [14(12)] of finite multiplicity which can be represented by means of the 


Min-Max Principle: 


An[E(I2)] 


inf sup 
EcV{n) ueE 
dim E=n 


Qnju) 


Correspondingly, there exists an orthonormal basis in L^(I2) of eigenfunctions [14(12)] 
associated with the eigenvalues A^[14(12)]. 

Note that, since the coefficients are fixed and bounded then 


QK^iu) < C'||u||^m,2(f^), 

for all u G II4”^’^(I2) where C is a positive constant independent of u and 12. Thus, 

1 /2 

since we assume that the space (14(12), (•)) is complete, we have that 


C||u||v[/m,2(f7) < QII^{ 


for all u G 14(12), where c is a positive constant independent of u. In other words, the 
two norms (•) and jj • equivalent in E(I2). Note that in general the 

constant c may depend on 12. However, stronger assumptions on the coefficients allow 
us to get c independent of 12. For example, if the coefficients satisfy the uniform 
ellipticity condition 

|a| = |^|=m |a|=m 

for all X G ^ = (CQ!)|a|=m G then it is straightforward that c can be chosen 
c = min{\/0,1} which is independent of 12. Condition (j2.7p will not be used in Section 
2 which is devoted to a general stability theorem. However, we shall use it in the 
following sections devoted to applications. 
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2.2 £^-compact convergence 

Let 0 be a fixed open set and V (fi) its corresponding space as in the previous section. 
For all sufficiently small e > 0 we consider perturbations fdg of and we denote by V (n^) 
the corresponding spaces of functions defined on We assume that the coefficients 
Aa/s are fixed functions defined in the whole of and that the operators and 

are well-defined and have compact resolvents. 

We denote by £ the extension-by-zero operator, which means that given a real¬ 
valued function u defined on some set in £u is the function extended by zero 
outside the given set. Clearly, for each e > 0, £ can be thought as an operator acting 
from L^(n) to consisting in extending the function by zero to all of and 

then restricting it to fig. As a matter of fact, this operator £ will be the key to compare 
functions and operators defined in fd and fig. The following concepts and definitions go 
back to the works of F. Stummel, see [33] and G. Vainniko see [34[ |35| among others. 
We also refer to Ham. Here we denote hy C{X) the space of bounded linear operators 
acting from a normed space X to itself. 

Definition 2.8 i) We say thatv^ G L‘^{£le) £-converges to v ^ LA{Q) if\\Vf: — £v\\i 2 ^Q^^-^ 
— >■ 0 as e — )■ 0. We write this as v. 

ii) The family of bounded linear operators G C{LA{£1^)) ££- converges to B ^ 
C{L‘^{£1)) if B^Ve Bv whenever v^ v. We write this as B^ ^ B. 

in) The family of bounded linear and compact operators B^ G C{L'^{Qe)) £-compact 

ss 

converges to B G £(L^(f7)) if B^ ^ B and for any family of functions v^ G L?‘{£Lf:) 
with < 1 there exists a subsequence, denoted by v^ again, and a function 

w G L?‘{£l) such that B^v^ w. We write B^ B. 


There is a strong relation between the T-compact convergence of a family of oper¬ 
ators and their spectral convergence. By this, we mean the convergence of eigenvalues 
and the associated spectral projections, see O Section 2.1]. Since in this particular 
work we are mainly dealing with B and Bg which are the inverses of the operators 
and defined above, we will define the spectral convergence just for this 

special type of operators. Hence, if we denote by the eigenvalues and 

eigenfunctions of and by {(A^, (/>n)}^i the eigenvalues and eigenfunctions of 

, (where we understand that the eigenfunctions are extended by zero outside 
and £l and they are normalized in L^(]R'^), we will say that the spectra behaves con¬ 
tinuously at e = 0, if for fixed n G N we have that A^ —^ A^ as e —)• 0 and the spectral 
projections converge in that is, if a 0 {An}^i, and Xn < a < A„+i, then if we 

define the projections —)• then 

sup{||P^(V’) — -Pa (V’)IIl 2 (]rjv) : fj G ||V’IIl2(rJV) = 1} 0, as e — )• 0. 

The convergence of the spectral projections is equivalent to the following: for each 
sequence 0 there exists a subsequence, that we denote again by e^, and a complete 
system of orthonormal eigenfunctions of the limiting problem {f)n})fLi such that — 
'Pn\\L^{MN) —>■ 0 as /c —>■ OO. 

As a matter of fact, we can show the following (see [U Thm. 4.1], |15]). 
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Proposition 2.9 Assume the operators satisfies the condition —)> ||it||L 2 (f 7 ) 

for each u € Lfi{Sl). If ^ £ C{L‘^{SLf)) are compact and ^ H~^^y then we 

have the spectral convergence of to Hyf^^y 

2.3 Examples 

We consider in this section some relevant examples of higher order operators. 

2.3.1 Polyhamonic operators 

An important class of higher order operators is given by the polyharmonic operators 
which we brefly discuss here as a prototypical example, see for instance [23] . 

For m € N, we set = baprnXja}. for all a,fi £ with |a| = |/3| = m, where 

^a/3 = 1 if a = /S and 6ap = 0 otherwise. With this choice, condition (12.7|) is satisfied. 

Let k £ No, 0 < A: < m and E(n) = n endowed with the norm 

(j2.ip of li k = m and Si has finite Lebesgue measure then E(n) is a closed 

subspace of W^’'^{Sl) and the embedding E(fd) Lfi{Sl) is compcict. If 0 
m then, under very weak regularity assumptions on (say, has a quasi-resolved 
boundarjtl in the sense of Burenkov |lll §4.3]), is a closed subspace of W'^’‘^{Sl) 

and V(fl) “—)• Wq if in addition fl has hnite Lebesgue measure then the embedding 
Wq{SI) LP‘{SI) is compact, hence the embedding E(n) Lp‘{Sl) is also compact. 
If A: = 0 and the open set Si is bounded and of class see Definition ED below, then 
the embedding I/(D) ^ L?‘{Sl) is compact (see Burenkov HD Thm. 8, p.l69]). 

Note that if A: = m then E(D) = W^'‘^{Sl) and integrating by parts one can realize 
that 

A~uA~vdx, if m is even, 

VA^“ uV A^~ vdx, if m is odd, 

for all u, u £ W^'‘^{Sl). In this case we obtain in (12.6|) the operator, Lu = (—A)™'u + u 
subject to the Dirichlet boundary conditions u = ^ = ... = 0 on dSl. Here 

and in the sequel u denotes the unit outer normal to dSl. The operator (—A)™'?/ is the 
classical polyharmonic operator of order 2m. 

In the general case k < m, the classical problem reads 

{—A)^u + u = f, in D, 

gL = 0, V j = 0, ...,A;-1, ondn, 

BjU = 0, V j = 1,..., m — /c, on dSl, 

where Bj are uniquely defined ‘complementing’ boundary operators. See Necas [31] for 
details. 


E 

|o| = |/3|=m 


Aai 3 D°‘uD^v dx 


^this includes the case of open sets satisfying the classical cone condition, as well as the case of open 
sets of class C°, see Definition ITT] below. 
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2.3.2 Biharmonic operator with Dirichlet, Neumann and Intermediate bound¬ 
ary conditions 


Let us consider the case m = 2 in the previous example. The quadratic form (12.31) will 
read 


/ D'^u : D^vdx + / uvdx, (2-10) 

Jn Jn 

for all u,v in an appropriate energy space lL(n). Here and in the sequel D^u denotes 
the Hessian matrix of u and D^u : D^v = Yl^j=i dxidx ■ dxidx ■ Frobenius product 

of the two matrices. 

As above, if V(fl) = bFgthe classical operator (j2.5p is given by the biharmonic 
operator and we obtain the classical Dirichlet problem 


A^u + u = f, 

u = 0, 

^ — n 


in Q, 
on dQ, 
on dQ. 


It is well-known that if N = 2 the Dirichlet problem for the biharmonic operator is 
related for example to the study of the bending of clamped plates. 

If V(fi) = IF^’^(H), by using the ‘Biharmonic Green Formula’ (I8.58P we obtain the 
classical Neumann problem 


{ A'^u + u = f, in H, 

0 = 0 , on on, 

divgQ((D\) ■ i')gQ -1-^ = 0, on dD, 

involving the well-known tangential divergence operator, see Section f8.4l for basic def¬ 
initions. We recall that if N = 2 the Neumann problem for the biharmonic operator 
arises for example in the study of the bending of free plates. See also Chasman |18] . 

Finally, If H(D) = n IFq^’^(H), proceeding as above we obtain the classical 

intermediate problem 


A^n + u = f, in H, 

u = 0, on dQ, (2.11) 

0 = 0, on dQ. 


We recall that if A = 2 the intermediate problem for the biharmonic operator arises 
for example in the study of the bending of hinged plates (sometimes called simply- 
supported) . 

We note that since u = 0 on dfl then the second boundary condition in (j2.1ip can 
be written as 


Au-aA=0, 

ov 

where K is the mean curvature of the boundary, i.e., the sum of the principal curvatures. 
See Gazzola, Grunau and Sweers [Ml for further details. 


3 A general stability theorem 

The following condition on and H(He) will guarantee that Hy^Q ^ is T-compact con¬ 
vergent to in the sense of Definition 12.81 
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Definition 3.1 (Condition C) Given open sets ilg, e > 0 and 12 in and cor¬ 
responding elliptic operators defined on respectively, we say that 

condition (C) is satisfied if for each e > 0 there exists an open set ilg C 12 H 12^ such 

thaE 

lim 112 \ KJ = 0, (3-2) 

£->■0 

and the following conditions are satisfied: 


(Cl) If v^ € l^(12e) and sup^yQQQfiv^) < oo then 

limJveW L^n,\K,) = 0; (3-3) 

(C2) For each e > 0 there exists an operator : 1^(12) —)• 1^(12£) such that for all 
fixed if £V (12) 


(i) limQKfiT.ip-if) = 0, 

(ii) = 0, 
£—>■0 

(iii) sup||r,(^||i 2 (f^^) < oo. 
£>0 


(C3) For each e > 0 there exists an operator E^: from y(12e) to such that the 

set £'e(i/(12e)) is compactly embedded in L^(12) and such that 


(i) If Ve G 1^(12£) is such that sup^^g Qfie(^£) < oo then lim QxfiE^v^ — vf) = 0; 

£—>■0 


(ii) sup sup 


\EfV 




>o„ey(C,) QJ (v) 


< oo; 




(iii) // Ue G l^(12e) is such that sup^^g Qr 2 e('(^£) < oo and there exists v G L^(12) 
such that, possibly passing to a subsequence, we have \\E,,v^ — t'||L 2 (Q) 0, then 

V G i/(12). 


Example 3.4 Consider the simpler case 12 C 12e and 1/(12) = l/(12e) = 

ll/”*’^(12e) for all e > 0, and assume Q is satisfied. We set ilg = 12 for all e > 0. 
Thus II.SM is trivially satisfied. Assume that 12 is sufficiently regular to guarantee the 
existence of a bounded extension operator from W'^’‘^{Tl) to 11/”*’^(M'^). Then condition 
(C2) is satisfied: indeed, the extension operator may serve as operator T^. As far as 
the operator E^ is concerned, one can use the restriction operator: in this way, the 
compactness of the embedding W'^’‘^{TL) ^ lfi{TL) allow to conclude that condition 
(C3) is satisfied . Thus, in order to verify the validity of condition (C) it just suffices 
to check the validity of (Cl). 


We now prove the following general statement. 


Theorem 3.5 Let £ be the extension-by-zero operator. If condition (C) is satisfied 
then 4 H-\ 

^this condition guarantees that || 2 ^u||L 2 (n^) —>■ ||w|L 2 (n) which is the basic hypothesis in the theory 
of S-convergence and it is an assumption of Proposition 12.91 





Proof. First we prove that ^ Let be a sequence £i-convergent 

to / G L^(r2), i.e., 

VimJf,-£f\\L2^n,)=0. (3.6) 

Let Ue G F(r2e) be such that 

Hn.v, = /„ (3.7) 

and let u G V{Q) be such that Hqu = /. We have to prove that Ue is £i-convergent to 
u, i.e., 

lim ||ue^ 0. (3.8) 

We will prove this statement by showing that for any sequence —>• 0 there is a 
subsequence e(, —>• 0 which satisfies (13.8p . Moreover, in order to avoid a complicated 
notation with too many indices and subindices, we will keep denoting the sequences 
and subsequences by e. 

By (j3.7p and the Holder inequality it follows that 

QnAve) = {fe,Ve)L^n,) < IIIIl2 (o,)QJ//(^^ e)• (3-9) 

Since is ^’-convergent to /, there exists M > 0 such that ||/e||L 2 (r 2 £) ^ ^ for ah 

e > 0. Thus by (j3.9p it follows that 

supQo^(ue) < oo. (3.10) 

£>0 

By condition (C3) (ii), it follows that is a bounded sequence in iy"*’^(H). Accord¬ 
ingly, by the compactness of the embedding E^(V{£l^)) C L?‘{£t) and the reflexivity of 
the space W”^’^(H) there exists u G VF™'’^(H) such that, possibly considering a subse¬ 
quence, Ef^Vf: converges to u strongly in and weakly in iy”*’^(H) as e —)• 0. By 

(13.101) and condition (C3) (iii) it follows that u G H(H). We now prove that u = u. Let 
G H(H) be fixed. Since the operator takes values in H(He) we can use T^ip as a 
test function in the weak formulation of the problem in He and obtain 

QnAve,Teip) = (3.11) 

It is easily seen that 

QnA'>Je,Te(p) = + Qn,\KS'^e,T^^) 

= QrAv^,^) + QKAve,Teip - (f) + 

= QxAEeVe,^) + “ E^V^,^) + “ 7 ^) + Qn,\KA'^e,Teip) 

= QQ,{E,,Ve,(p) - QQ\KS^eVe,^) + Qk, (t’e “ E^V^,^) + QKAVe,T^^P “ 
+QOe\X£(r’e)Le7^)- (3-12) 

By the boundedness of the coefficients Aa/s, the space (^”^’^(12), || • \\w^^‘^{n)) is 
continuously embedded in {Q), {■)). It follows that the sequence HgUe is 

weakly convergent to u in {Q), {■)), hence 

lim Qn{E^Ve, 7 ?) = Qq(u, ^). (3.13) 

e-i^O 

Since QQ\KSEtVf:) is a bounded sequence, if is fixed and lime_j.o |H \ = 0 it is 

easily seen that 

limQn\X,(-EeUe,V7) = 0. (3.14) 
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Since QkX^) is a bounded sequence, by condition (C3) (i) it follows that 


lim QK.i'i’e - E^Ve,^) = 0. 


e^O 


(3.15) 


By (13.101) Qxei'^e) is a bounded sequence hence by (C2) (i) it follows that 


Similarly, by (C2) (ii) 


-(/?) = 0. 

e^O 


= 0 . 

e-)-0 ^ 


Thus, by (I3.12|) - (I3.16D it follows that 

lim Qq^{v^,T^ip) = Qniu^if) 
e^O 

Moreover, 

(/ej = (/e — £ f ,T^ip) + {£ f ,T^ip) ^2 

By (j3.6p and condition (C2) (iii) we have 

lim(/e -T/,reV?)i2(Q^) = 0. 


(3.16) 

(3.17) 

(3.18) 

(3.19) 


Furthermore, 


l(^/)T’£99)2,2(n^) - (/, (/?)2,2(Q)1 < \{f,T^ip - ip)L2(^K^)\ 

+ \if L2yQ^Q^)\X^) \ + \if,^)L^{n\Ke)\- (3.20) 

By condition (C2) (i) the first summond in the right-hand side of (13.201) vanishes as 
e —)■ 0. Moreover, since |(n H fl^) \ Kf:\ —>• 0 as e —)■ 0, by condition (C2) (iii) the second 
and the third summonds in the right-hand side of (13.201) vanish as e —>• 0. Thus 

lim(T/,T,(^)i2(f2^) = if,ip)L2{Q)- (3.21) 

By (j3.1ip . (j3.17p . (j3.18p . (j3.2ip it follows that 

Qniu,ip) = {f,ip)L 2 (^f^y yipeV{n). (3.22) 

Since, u € 17(fl), we have u = u. 

Observe now that 


= Il^elli2(n,\n) + ll^e “ '“llL2((n,no)\i<',) + ll'^e “ ^IIl2(^^).(3.23) 

Note that 


II^C — \\EeVe ^Ili2(/f^) + \\EeVe 'l’e|li2(j:4'^) 

< \\EeVe — u\\j^2^yi'^ + \\Ef^Ve — (3.24) 

Moreover, 

II''^c|Il2(q^\ 0) < ll^e|li2(Q^\;^^), (3.25) 
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and 


ll^c ^llL2((n^nn)\A'd — + \\''J'\\‘L^{n\Ke)- (3.26) 

Thus, by (Cl), (C3) (i) and (I3.23ll - (l3.26p it follows that (13.81) holds. Thus, is 
T-convergent to as e ^ 0. 

Exactly the same argument can be used to prove that if /e G L^(fle) is such that 

SUp||A||L2(f^^) < OO 
£>0 


and De G E(ne) is such that then there exist u G L^(n) such that, possibly 

considering a subsequence. 


lim llue 
£^•0 


^^IIl2(q^) — 0. 


(3.27) 


This implies that is £^-compact convergent to ^ as e ^ 0. 


□ 


4 Dirichlet boundary conditions and Mosco convergence 

In this section we consider the operator (12.61) on a bounded open set 12 in , subject 
to Dirichlet boundary conditions 



dm 1 ^ 


on (912. 


(4.1) 


This has to be understood in the general frame discussed in Section 2 as follows. 

Imposing Dirichlet boundary conditions to the operator L on 12 means that the 
domain V (12) of the corresponding quadratic form is given by 

V(12) = Wo™’^(12). 


This will be understood throughout this section. Here we assume that the coefficients 
Aap satisfy the uniform ellipticity condition (12.71) . We recall that if 12 is bounded then 
the Sobolev space VE™’^(12) is compactly embedded in L^(12). Thus, as it is explained 
in Section 2, the operator is well-defined and has compact resolvent. 

The spectral stability of higher order operators subject to Dirichlet boundary con¬ 
ditions on variable domains was discussed in Babuska and Vyborny [8] where sufficient 
conditions ensuring stability were given. Those conditions are nowadays understood in 
the frame of the notion of Mosco convergence which we now recall. 


Definition 4.2 Let D be a bounded open set in . Let e > 0 be a family of 
open sets contained in D. Let m G N and 12 be an open set in D. We say that the 
spaces converge in the sense of Mosco to the space as e ^ 0 if the 

following two conditions are satisfied: 

(Ml) For any ip G IT™’^(12) and e > 0 there exists ip^ G such that ip^ ^ ip in 

as e^O. 

(M2) If V £ W^’‘^{D) and there exists a sequence G W()”’^(12c^) such that v 

in W^’‘^{D) as —)• 0, then v G 
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Note that in the previous definition it is understood that functions ip, are ex¬ 

tended by zero outside their domain of definition. Moreover, the condition v G 
in (M2) has to be understood in the sense that the function v can be approximated in 
^ sequence of C°°-functions with compact support in 0. 

We prove the following expected result. 

Theorem 4.3 The Masco convergence of the spaces to as e ^ 0 

implies the validity of condition (C), hence the £-compact convergence of the operators 
the operator 

Proof. Assume that D, and 12 are as in Definition 14.21 We set ifTe = Dg n 12. 
We divide the proof in several steps. 

Step 1 We prove that condition (13.21) hold. Using standard properties of the 
Lebesgue measure, to prove (13.2p it is enough to show that for any compact set if C 12, 
we have [AT \ 12^1 —)• 0. But, since K is compact and K C 12, we have the existence of 
a function p G C“(12) with p = 1 m. K. From (Ml) we have a sequence of functions 
Pe G VF™’^(12e), such that p^ ^ p m. and in particular in L^{D). But this 

implies that 

|{x G 12 : \peix) - p{x)\ > 1/2}| 0 

and therefore, since in A' \ 12^ we have p(x) — Pe{x) = 1, then A' \ C {x G 12 : 
\pe[x) — p{x)\ > 1/2} and therefore \K \ 12^1 —)• 0. 

Step 2 We prove that condition (Cl) is satisfied. Assume that G 
are as in (Cl) and assume directly that are extended by zero outside Dg. Since 
Ve G VF™’^(D), by the Sobolev’s Embedding Theorem G L'P{D) for some p > 2. 
Thus, 

\\Ve\\L^{n\K^) < \^\ Ke\'^ P || ||Ue||< c| 12 \ AT^[ 2 p 

where c > 0 is a constant independent of e. This, combined with (13.21) implies that 
IIu||l2(o\K£) —)• 0 as e —)• 0. Thus, in order to prove that condition (Cl) is satisfied, it 
suffices to prove that ||i'e||L2(OE\o) ^ 0 as e —>■ 0. Assume by contradiction that this is 
not the case. Then there exists a subsequence v^^ such that ||'ye„ ^ c > 0. 

Moreover, possibly passing to a subsequence, by (M2) there exists v G 1U™’^(12) such 
that converges to v in L^{D). In particular, it follows that ||r'|| 2 , 2 (£)\Q) = c which 
contradicts the fact that v G W™’^(12). 

Step 3 We prove that the validity of (Ml) implies the validity of (C2). For any 
p G B/)"’^(12), we set T^p = pe where pe is as in (Ml). Obviously, we have that 

Q K,:{TeP ~ p) c||(/9e — p\\-^m,2^j^y 

where c > 0 is independent of e. Since p^^ ^ pin W^’‘^{D), follows that condition (C2) 
(i) is satisfied. We now prove that condition (C2) (ii) is satisfied. We note that there 
exists a constant c > 0 independent of e such that 

Qn^\KS^<i) ^ < c\\p^ - p\\w”^^'^{np\K^) 

+c\\p\\w^,'2{n^\K^) = \\Te - + h\\p\\wrn,2{p,\K^)- ( 4 - 4 ) 

It is now clear that by (Ml) and the absolute continuity of Lebesgue integrals, the 
right-hand side of (14.41) goes to zero as e —)• 0, hence condition (C2) (ii) is satisfied. 
Condition (C2) (iii) is trivial. 
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Step 4 We prove that the validity of (M2) implies the validity of (C3). We set 
Ef:{v) = Extof where Exto is the extension-by-zero operator. Clearly, conditions (C3) 
(i), (ii) are trivially satisfied. We now consider condition (C3) (hi). Let be as 
in (C2) (hi). Then Exto^e is a bounded sequence in hE™’^(L)). Thus, there exists 
V G W^’ (D) such that, possibly passing to a subsequence, Extou^ is convergent to v 
strongly in L‘^{D) and weakly in VE™’^(Z)). By condition (M2), we immediately have 
that V G hence condition (C2) (ih) is also satisfied. 

Remark 4.5 Sufficient conditions ensuring the Mosco convergence of spaces 
to are well-known. We refer to Bucur and Buttazzo and Henrot 121^ for 

a detailed discussion in the case m = 1. We note that such conditions typically involve 
geometric notions describing the vicinity of sets (for example, the Hausdorff distance) 
as well as uniform regularity or topological assumptions on the domains. Some of the 
conditions known in the case m = 1 easily extends to the case m > 1, as in the case 
of the compact convergence of sets. For example, if C Cl is a sequence of open sets 
compact convergent to Cl as e ^ 0 (i.e., for any compact set K dCl there exists ck > 0 
such that K C Cl,, for all 0 < e < ex) then one can prove that the spaces W0^''^{Cl,) 
converge in the sense of Mosco to fL))"’^(n) as e —>■ 0. See Babuska and Vyborny 
for more information. 

5 Neumann boundary conditions 

In this section we consider the operator (j2.6l) subject to Neumann boundary conditions 
on bounded open sets hi in This has to be understood in the general frame discussed 
in Section [2] as follows: by Neumann boundary conditions we mean that the domain 
V (12) of the corresponding quadratic form Qq is given by 

V{Cl) = W^’\Cl), 

and this will be understood throughout this section. Here we assume that the coeffi¬ 
cients AajS satisfy the uniform ellipticity condition (|2.7I) . 

It is well known that both the smoothness of the domains and the kind of per¬ 
turbations that we are allowed when dealing with operators with Neumann boundary 
conditions is more restrictive than in the Dirichlet case. An appropriate setting for this 
issue is clarified with the notion of atlas, as for instance in m Definition 2.4]. Eor the 
sake of completeness and clarity, let us include here the definition. 

Eor any given set V G and h > 0 we denote by Vs the set {x G : d{x, dCl) > 
h}. Moreover, by a cuboid we mean any rotation of a rectangular parallelepiped in . 

Definition 5.1 [Definition 2.4, fWj] Let p > 0, G N with s' < s. Let also 
be a family of bounded open euboids and be a family of rotations in 

M.^. We say that A = (p, s, s', is an atlas in with parameters 

p, s, s', {Vj}j^i, {rj}j^i, briefly an atlas in Moreover, we eonsider the family of 
all open sets Cl C satisfying the following: 

i) Cl C and {Vflp n D / 0 

ii) Vj n dCl A 0 for j = 1,... , s' and Vj n dCl = 0 for s' < j < s 

Hi) for j = 1,..., s we have 

D'(^) = {a: G : atj < Xi < bij,i = 1,... ,N}, j = 1,... ,s 
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rj{Vj n S7) = {x G : a^j < xjv < 9 j{x)}, j = 1,... ,s' 

where x = (xi,..., X 7 v_i), Wj = {x G : a^- < Xj < bij,i = 1,..., — 1} and the 

functions gj G C^’'^(Wj) for j = 1,... ,s', with A: G N U {0} and 0 < 7 < 1. Moreover, 
for j = 1,..., s' we have onj + P < djix) < b^j — P, for all x G Wj. 

We say that an open set Q is of class if all the functions gj, j = 1,..., s' 

defined above are of class C^’'^{Wj) and \\gj\\c’^,'i(Wj) — open set 

O is of class if it is of class {A) for some M > 0. Also, we say that an 

open set Q is of class if it is of class {A) for some atlas A and some M > 0. 
Finally, we denote by the class for fc G N U {0}. 

We recall that if fl is a bounded open set then the Sobolev space is 

compactly embedded in see e.g., Burenkov [TT]. Thus, as it is explained in 

Section 2, the operator is well-defined and has compact resolvent. 

In this section we discuss the T-compact convergence of the operators ) 

on families of open sets e > 0. The following theorem is in fact a generalization to 
higher order operators of [T], [3l Prop. 2.3]. 

Theorem 5.2 Let Q be a bounded open set in M.^ of class and Hg, with e > 0, be 
bounded open sets in of class C^. Let A be an atlas in , as in Definition 15.11 
and M > 0. Assume that for each e > 0 there exists an open set K,. C LI Ci Ll,, of class 
C^^{A) satisfying 113.^) and one of the following two equivalent conditions: 

i) If Ve G and sup < 00 then lirn = 0; 

ii) limTe = 00 , where 


Then we have limjric \ = 0 and condition (C) is satisfied. Hence, H„}^ 2 (r, \ ^ 

£—>■0 

ij—i 

Proof. First we note that conditions i) and ii) both imply that lim jlle \ iPej = 0. 

€—^0 

The proof of this is very similar as the one from [3l Prop. 2.3] and we skip it. 

That i) implies ii) is very simple. Indeed, if is bounded for some sequence —)• 0, 
then one can find functions G W'^’‘^{Llek) with \\veJ\L2{n^^\K^i^) = 1 and Qn.^ivek) 
bounded. Hence, jjuej, ) is a bounded sequence, which contradicts i). 

We now prove that ii) implies i). If i) does not hold, one can find a sequence 0 
and functions G with < 00 and = 

1. Since is of class C^{A) for all e > 0, there exists a bounded linear extension 
operator Exti^^ : —)• IT'”’^(M'^) with a uniformly bounded norm, i.e., 

SUp||Ext/C£ ||^^m,2/^^\_^^m,2mJV\ <00; (^-3) 

see Burenkov mi Thm. 3, Chp. 6] or the classical Stein’s book |36| . 


Tf = 


inf 

„gVV'"'2(Q,)\{0} 
11=0 on Ke 


QnAx) 
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We set 14^ = Extxj (I’efe )• By using the Sobolev’s Embedding Theorem and 
(1531) one can prove as in [3] that 0- Consider now the func¬ 
tion Wet, = - Veu- B is clear that w^t, = 0 on Ke^,, sup^^N (^e/c) < oo and 

liminffc_j,oo ll'u^ej. > 0. It follows that sup^gpjTej, < oo which contradicts ii). 

It remains to prove that condition i) implies that condition (C) is satished. By (12.71) 
it follows that the norm is equivalent to the Sobolev norm (12.11) . Thus condition 
(i) implies the validity of condition (Cl). 

Since 11 is of class there exists a bounded linear extension operator Exto : 
IT”^’^(n) —)• iy”*’^(M'^). Eor every e > 0, let be the operator of to 

dehned by = (Extn(/9)|Q^, for all (p G Since = p on and T^p = 

ExtQ<y9 on He for all p G IT"*’^(n), and lime_).o |lle \ Ke\ = 0, it follows that satisfies 
condition (C2). 

For every e > 0, let be the operator of to iy”*’^(n) dehned by E^u = 

Exti^^(tt|j^J, for all u G It is obvious that condition (C3) (i) is satished. 

Moreover, by (|5.3I) also condition (C3) (ii) is satished. We now prove that (C3) (iii) 
is satishes as well. Let G IT™'’^(lle) and v G L^(ll) be as in (C3) (iii). Since 
is uniformly bounded, by (15.3p it follows that Exti^^(rt|j^J is uniformly 
bounded in VE™'’^(]R'^). By the rehexivity of 1T”^’^(M'^) there exists v G IT™'’^(M'^) 
such that Ve converges weekly to v in IT™’’^(M^) as e —>• 0. Thus, E^v^ converges weekly 
to V in L^(ll) as e —)• 0 hence v = v and v G IT™'’^(11). Thus, also condition (C3) (iii) 
is satished. The proof is complete. □ 

We have the following result, which considers the particular case where hi C 

Corollary 5.4 Let fl be a bounded open set in of elass and with e > 0, 
be bounded open sets in R'^ of elass with hi C hie. Assume that one of the following 
two equivalent conditions is satisfied 


i) Ifve G IT™-’2(hle) and sup,>o \\ue\\w^^^ {np < oo then lime^o ll?^jL2(OeW) = 0- 

ii) lim inf = oo. 

u=0 on ^ ^ 


Then we have hme_).o |hle \ hl| = 0 and condition (C) is satisfied. Hence, ) 

is £-compact convergent to • 

Proof. The proof immediately follows by setting = £l for all e > 0 and applying 
Theorem 15.21 □ 


We also can prove a simple criterion ensuring spectral stability for Neumann bound¬ 
ary conditions, which generalizes the condition given in Arrieta and Carvalho [3l § 5.1]. 
This criterion can be easily formulated in terms of the notion of the atlas distance 
which is introduced in m- 

Definition 5.5 (Atlas distance) Let A = (p, s, {14}j=i, {?’j}j=i) be an atlas in 
R^. For all hli,hl 2 G C'™'(A.) we set 

d^^(hli,hl 2 ) = max sup sup \D'^gij{x) — D'^g 2 j{x)\ , (5.6) 

j—l,...,s Q<|Q,|<m {x,X[,r)£rj(Vj) 
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where gij, g 2 j respectively, are the functions describing the boundaries of respec¬ 

tively, as in Definition \5.1[ 

Moreover, we set and we called ‘atlas distance’. 

The atlas distance clearly depends on the atlas but has the advantage of being 
easily computable. In the case of open sets of class C^, dj^ is equivalent to the usual 
Hausdorff distance. Amongst the basic properties of it is also worth mentioning 
that is a complete metric space. See [13] for more information. 

Then we can prove the following 

Theorem 5.7 Let A be an atlas in and Q be a bounded open set of class C^^^{A). 
Let rie, e > 0, 6e bounded open sets of class C{A) such that 

lim dji,{Llf^, Ll) = 0. 


Then condition (C) is satisfied, hence ^ 

Proof. Let iLg be bounded open sets of class with M > 0 independent 

of e, such that C 0, Cl and such that < 2d_A{Lt,,,,LL). We denote by 

j and gKe,j the functions describing the boundaries of fie and respectively, as in 
Definition 15.11 

Recall now that if a function / belongs to a Sobolev space the type lT™'’^(a, 6) 
where (a, h) is a bounded real interval and /^^(a) = 0 for any i = 0,..., m — 1 then the 
following Poincare inequality holds 

ll/llLRa,6)<C’(6-an|/(™)||i2(,,,), (5.8) 

where C > 0 depends only on m. 

Let V G lP™'’^(De) be such that = 0. Then by applying Fubini-Tonelli’s Theo¬ 
rem, using (15.81) and the notation from Definition 15.11 we have 


\L^{n,\K,) ^ 


til 




=1 JWj JgKe.ji^) 


\v{x,X]s[)\‘^dxNdx 


p r 

< \9n,,j{x)-gK,,j{xl'^ 

j=l JWj Ja 

< Cd^ (Kf:,Ll,,^)\\v\\^rm,2|^Qy, 




dmy 


dx 


N 


dxridx 


(5.9) 


by which we immediately deduce that hme_>.oTe = oo, where is defined in Theo¬ 
rem 15.21 Thus, Theorem 15.21 allows to conclude the proof. □ 


Remark 5.10 ITe note that in Theorem 5.7 the assumptions on the open sets are 
quite weak. Indeed, it is not required that the sets belong to a uniform Lipschitz 
class and it is only required that they are of class C (.4.). In particular, the modulus of 
continuity of the functions describing their boundaries may blow up as e —)• 0. 
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6 Intermediate boundary conditions 

In this section we consider the operator (|2.6I) subject to intermediate boundary con¬ 
ditions on bounded open sets H in with smooth boundaries. By intermediate 
boundary conditions we mean that the domain of the corresponding quadratic 

form Qo is given by 

v{n) = n Wo’^{n), (e.i) 

with m > 2, where A: G N, 1 < A: < m, is fixed. This will be understood throughout 
this section. By well-known estimates for intermediate derivatives (see Burenkov m 
p. 160]) it follows that is a closed subspace of IT™'’^(n). 

Here we assume that the coefficients Aa/s are hxed and satisfy the uniform ellipticity 
condition (|2.7p . Thus the operator TAym) is well-dehned and has compact resolvent 
since H(n) is compactly embedded in Lr{Q). 

In this section we discuss the T-convergence of the operator ^ on suitable 

families of smooth open sets e > 0. Our analysis includes the case of open sets 
with oscillating boundaries. 

We consider, as in Definition 15.1] a fixed atlas A and constant M > 0 and assume 
D G CJ^{A). We will also consider that G C^^{A) for some constants not 
necessarily uniformly bounded in e. To simplify the proofs of the results, we will 
consider that the perturbation of the boundary is localized in just one of the cuboids 
Vi for some i = 1,... ,s', that is, one of the cuboids which touch the boundary. We 
refer to Corollary 16.181 for a general statement when the perturbation acts not in just 
a single cuboid. 

Therefore, let us denote the cuboid by V (we will drop the subindex) and hence 
we will assume that D \ V), = Dg \ V),, that is, the perturbation is localized in the 
interior of V. Without loss of generality, we may assume that V = W x {a,b), where 
IT = {x G : aj < Xj < bj,j = 1,... ,A^ — 1}. Moreover, T n D = {{x,xn) G 

ITx]a, b[: a < xn < gix)} and T n = {(x, xn) G ITx]a, b[: a < xn < 9eix)}. The 
functions g and define the boundary of D and in V and as in Definition 15.11 we 
assume that a+p < g,g^ < b-p, g,ge G C'”'(IT) with \\g\\c^(w) < \\ge\\c’^{w) < 

The following lemma provides a sufficient condition for the T-compact convergence 
of ^ to Hy^^y As usual, || • ||oo denotes the L°°-norm. 

Lemma 6.2 With the notation above and assuming that for every e > 0 there exists 
Kg > 0 such that 

(i) Kg > IlfiTg - V e > 0 

(a) limg_>.o Kg = 0 

(Hi) limg_j.o = 0, V /3 G with |/3| < m. 

kT ^ 

Then condition (C) is satisfied, hence ^ 

Proof. Our argument is based on the construction of a suitable diffeomorphism 
from Dg onto D which coincides with the identity outside Vp. Thus, since Dg\I/), = D\I/), 
for all e > 0, we can carry out our construction in V and assume directly that D = DnT 
and Dg = Dg n V. Hence, D = {(x,XAr) G ITx]a, 5[ : a < xn < g{x)} and Dg = 
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{(x,XAr) S Wx]a,b[ : a < < 5 e(x)} where W is defined above. LefH k > 2{m + 1 ). 

We set fcg = and ge = Qe — ke, = {(xjXat) G Wx]a, 6[ : o < xn < 5e(®)}- Note 
that by (i), (ii), 0 7 ^ C 0 H for all e > 0 sufficiently small. Moreover, by (ii) 
condition (13.21) is satisfied. 

Throughout the proof, we will denote by C a generic constant which will be inde¬ 
pendent of e and all the functions involved. This constant may change from line to line. 
If at some point we want to distinguish some constant we will use another notation and 
make this clear. 

We now prove that condition (Cl) is satisfied. Let e > 0, be as in condition (Cl). 
We will need the one dimensional embedding estimates \\f\\L<^{a,b) < ^\\f\\w^’^{a,b) 
where the constant K = K{d) is uniformly bounded for \b — a\ > d (see, e.g., Bu¬ 
renkov m)- Hence, by Tonelli’s Theorem, and applying this last estimate to the 
function Ve{x, •) in the interval (a, 5 e(x)), we get 

/ \Vf:\dx= / / \Vf:{x,XN)\dXndx 

Jw Jgd^) 


< C 


/ \ge{x)-ge{x) 
Jw 


WAX,-] 




^\\dx 


(6.3) 


< <^1156 - ge\\L--(W)\\Ve\\wh2(nA ^ CkI QnAVe) ^ 0 . 

By (j6.3p the validity of condition (Cl) follows. 

We now prove that condition (C2) is satisfied. Let : Hg —)• H be the map dehned 
by $e(x,XAr) = {x,X]\f — h^{x,XN)) for all {x,xn) G where 


he{x,XN) 


0, CL < Xn < ge{x), 

{g,ix) - g{x)) , if ge{x) < Xn < ge{x) • 


(6.4) 


Geometrically speaking, the map <l)e transforms the “vertical segment” 1% = {{x,xi\[) ■ 
a < Xn < gt{x)} to the segment = {{x,xn) ■ a < xn < g{x)}-, leaving invariant the 
part of the segment with a < xn < ge{x)- 

Note that for e fixed the transformation is a diffeomorphism of class C™ from 
He onto H. The C”^ norm of ‘be will not be bounded in general for m > 2 as e —)• 0, 
but by the choice of k and having in mind (i) (ii) and (hi), it follows that 


C-^ < |detT>$e| < C. (6.5) 

To see this, just note that |detZl<I>e| = |1 — dxffhe\. 

In order to estimate the derivatives (up to order m) of the transformation <I>e, we 
need to study the derivatives of the function h^- By the Leibniz formula we have 

XN - ge{x) ^ 


D°he{x) = 

0<7<a 


a 


7 


D^g,{x) - g{x))D'^-'^ 


ge{x) -h{x) 


for all {x,xn) G W x]a, b[ with g^{x) < xn < gt{x)- By standard calculus, it is easy to 
check that 

XN - ge{x) ^ 


jja--y 


geix) - geix) 


< 


c 


< 


c 


|g,(x)-5,(x)|l“l-H - >I-H’ 


( 6 . 6 ) 


®the value of k does not play any significant role and is used only to prove (1531) 
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hence, 


(6.7) 




< C 


E 

0<7<Q: 


ll^"(5e-. 


He 


|a|-|7| 


for all e > 0 sufficiently small. 

Let Tg be the map from V{^) to defined by 


Te(p = (po^^, 

for all (f G 17(0). Note that is well-defined since <l>e is a diffeomorphism of class C™. 

Condition (C2) (i) is immediately satisfied since = (p on for all p> G 17(0). 
We now prove that condition (C2) (ii) is satisfied. Let p) G 17(0). By the chain rule, 
for any multindex a with |a| = m, we have 

D“(^{#.(x))) = ^ Df‘^(^,{x))pZ,pi,) (6.8) 

l<|/3|<m 


where fii^e) is a homogeneous polynomial of degree |/3| in derivatives of of order 
not exceeding m — |/3| -|- 1, and coefficients depending on a but not depending on e. 
Since ‘^^(x) = x — (0, h^{x)), then a derivative of of order \(5\ is either constantly 1 or 
0 or a derivative or of order les or equal than |/3|. Then /^(‘be) is a polynomial of 
degree less or equal than |/3| in the derivatives of of order not exceeding m — |/3| -|- 1. 

In particular, using (j6.7jl this implies that 


IK,s(‘i>.)ll=o < c 


E 

0<|7|<m-|/3|-)-l 


ll^^(ge-g)||oo 

m-|^|+l-|7| 



l/3| 


Note that, in particular if |/3| = 1 then 




(6.9) 


, 0<|7|<m 


where, as it is customary, we denote by o(l) a function which goes to 0 as e —>• 0. 
If 2 < \j3\ < m, we have 


IK.,s(‘i>.)ll=o<c ^ 


11079. - 9)11. 


^0<|7|<m-l 

where we have used hypothesis (iii) 


m—1—I 7 I 


+ 1 <c 


( 6 . 10 ) 
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Hence, we have 


Qn.\KSTei^)) < f \H<^e)\^dx + C [ |Z)“((^(<h,))|2dx 


<C [ \<f\‘^dx + C Y [ \D^V’i‘^^ix))p^^p{<^e)\‘^dx 

J^\Ke Jn,\K, 

l<l/3|<m 



l<|/3|<m 


Notice that since (p € VH™’^(r2) is a fixed function and \ Xel —)■ 0, then 

^ ^ ( 6 - 12 ) 

Also, notice that since p G ^ Ty^’^(n), we have that for all i = 1,..., A^, 

Vp G Therefore, Vp{x,-) G lT^’^(a,^(a:)) L°°{a,g{x)) a.e. x G Vh. With 

a similar argument as the one we use to show (j6.3l) . we have 

ll¥^ll^l.2(r2\K£) — ^Wde — <?I|oo||aII^2 .2 (q) < Ck^. (6.13) 

Hence, from (j6.12j) we have that the first term in (jfi.lljl goes to 0. Moreover, for 
the second term in (j6.1ip we consider the sum for \I3\ = 1 and 2 < |,0| < m separated, 
apply (j6.9jl and (jfi.lOp to obtain 

y~! \\Pmp{^<d\\^cxi ll¥’ll^l/5|.2(Q\i^^) < o{1)k^ ^\\v\\w'^.^{Q\K^) + ^\\P\\w^<‘^{Q\K^) ^ 0 

\a.\=m 

l<\P\<m 

where we use (j6.12l) and (I6.13p . This shows that (C2) (ii) is satisfied. 

Condition (C2) (hi) is trivial. 

We now prove that condition (C3) is satisfied. By condition (iii) it follows that 
llV^ellcxD is uniformly bounded for e sufficiently small (recall that m > 2). Thus the 
open sets fig belong to the same class C\^{A) for a suitable fixed M > 0. Hence, there 
exists a bounded linear extension operator Ext^^ from to such 

that 

SUp||Exto^ < oo. (6.14) 

We set E^u = (Extf 2 ^n)|Q for all u G ^(0^). It is straightforward that satisfies con¬ 
ditions (C3) (i), (ii). We now prove that condition (C3) (iii) is satisfied. Let v^-, e > 0, 
and V be as in condition (C3) (iii). Since E^^v^ is bounded in W'^’‘^{Vl) and IT™'’^(fI) is 
compactly embedded in L^(n), it follows that there exists v G IT™'’^(I1) such that, by 
possibly considering a subsequence, converges weakly to v in IT™'’^(n) and strongly 
in L‘^{Q) as e —)• 0. It follows that v = v, hence v G IE™'’^(n). It remains to prove 
that V G Wq‘^{^). To do so it suffices use the extension-by-zero operator £. Indeed, 
£{v^) G Wq'‘^{U) where 17 is a bounded open set containing all sets and fl. Since 
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is bounded in Wq’‘^{U) and Wq’‘^{U) is compactly embedded in L‘^{U) it follows that 
there exists v G VKq’^(O) such that by possibly considering a subsequence converges 
weakly to v in Wq{U) and strongly in L'^{U) as e —>■ 0. Clearly, v vanishes outside Q, 
hence v G Obviously, D = u in 0. The proof is complete. □ 

Remark 6.15 It is not difficult to see that the fact that the set W is an {N — 1)- 
dimensional cuboid of the form IT = {x G : aj < xj < bj,j = 1,... ,N — 1} is not 

essential at all in the proof of Lemma \6.^ As a matter of fact, exactly the same proof 
works if we consider a general smooth, say piecewise , set W and V the cylinder of 
base W that is V = W x {a,b). 

Remark 6.16 By the classical Gagliardo-Nirenberg interpolation inequality ||T)^/||oo < 

C(I]|Q:|=m (cf e.g., ISE, P- 125]), it turns out that in order 

to verify condition (Hi) in Lemma WIE. it suffices to verify it for \ j3\ = 0 and 1/3| = m. 

We can deduce now the following. 

Proposition 6.17 With the notations above if Wg^ — ^Hoo 0 and if 
sup {\\ge - g\\^~^ WD^ge - D^^gWoo} ^ 0 , 

\a\=m 

then condition (C) is satisfied. Hence Hy^^ ^ ^v\n)- 

1 

Proof. Let us denote by (5e > 0 a sequence such that sup|„|=^{|| 5 e— \\D°‘g^ — 
D°‘g\\oo} < de 0 and define 

p, = max{(5e, \\g^ - g\\fff ~^} ^ 0 

Let us choose 

II5.-5II#" 

Ke = - 

Pe 

and assume directly that 7 ^ 0. Then, since 2/(2m — 1) < 1, \\g^ — p||oo —> 0 and 
Pe 0 , for e small enough we have 

\\ge - g\\oo < llffe - g\\^~^ < He< he “ g\\^~^ ^ 0 

and therefore hypothesis i) and ii) from Lemma 16.21 hold. 

Moreover, with the definition of k^, 6^ and and noting that Pe > 5e, we have 

sup - D^gWoo} < ph^'^de < d^^ ^ 0 

\a\=m 

Using now Remark 16.161 we can easily show that iii) from Lemma 16.21 holds and the 
proposition follows. □ 

In case the perturbation does not act only over one cuboid or it is different from 
the set V considered in Remark 16.151 above, we can also prove the following result with 
the aid of a partition of unity subordinated to the family of cuboids We will 

denote the functions that define the boundary of kl and fig in Vi for alH = 1 ,..., by 
Pi and g^^i respectively. 
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Corollary 6.18 With the notations above, if we assume that for every e > 0 there 
exists Ke > 0 such that 

(i) Ke>\\ge,i-gi\\cx>, V e > 0, Vz = l,...,s'. 

(a) lim£_>.0 Ke = 0 

(in) lim£_j,o = 0, V /3 £ with |/3| < m and for all i = 1,... , s'. 

Then condition (C) is satisfied, hence ^ 

Remark 6.19 Observe that a similar observation as in Remark \6.1(A can he applied in 
this case. 

Finally, we can deduce the following 

Theorem 6.20 Let A be an atlas in , M > 0, m £ N, m > 2. Let fig, Ll £ C(0{A), 
e > 0, be such that 

lim d^fi (fie, Q) = 0 . 

Then condition (C) is satisfied, hence ^v\n)' 

Proof. Apply Corollary 16.181 with = (d^^ . □ 

7 The biharmonic operator with intermediate boundary 
conditions 

In this section, we shall consider the biharmonic operator subject to intemediate bound¬ 
ary conditions in a family of domains with oscillatting boundaries. 

Without loss of generality and to simplify the exposition, let us assume that our 
domain kl C is of the form Lt = W x (—1,0) where W C is either an (iV — 1) 

dimensional cuboid as in Lemma 16.21 or a smooth domain as in Remark 16.151 We also 
assume that the perturbed domain is giyen by = {(x, xat) : x £ W, — 1 < xn < 
ge{x)} where 5 e(x) = e" 6 (x/e) for all x £ IF and b : —)• [ 0 , is a F-periodic 

hxed function, where Y is the unit cell Y = (— 5 , and a > 0. Note that for 

simplicity, we haye assumed that b > 0. Note that the general case can also be treated 
in a straightforward way, although to ayoid annoying technicalities and for the sake of 
the exposition we will stick to the simpler case b > 0. This implies in particular that 
C rie, e > 0. We denote by Fg the set {(x,XAr) : x £ IF, xn = ge{x)} which is the 
part of the boundary of aboye IF. It is also coyenient to set LIq = LI. In the sequel 
we shall also identify Fq = IF x {0} with IF. 

Namely, we shall consider the operators 

+ u 

on the open sets with u subject to the classical boundary conditions 

du 

u = 0 , and Au — K— = 0, on dLle, 
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where K denotes the mean curvature of i.e. the sum of the principal curvatures. 
More precisely, the operators are the operators associated with the quadratic 

form 

Qn^{u,v)= f {D^u : D'^v + uv)dx (7.1) 

defined for all n as discussed in Section 2. Recall that D^u 

denotes the Hessian matrix of u and D^u : D‘^v = 'Ylfj=i dTdx dx^dx • 

It is clear that this is a special case of those discussed in the previous section with 
m = 2 and A: = 1 in ()6.ip . 

It will be convenient to denote by the operator + / subject to Dirichlet 

boundary conditions on W, that is, 

u = ^ = 0, on W (7.2) 

ov 

and intermediate boundary conditions on \ W, which is the operator canonically 

associated with the quadratic form (17.ip defined for all u,u € IRq’^(H) where Wq’^{Q) 

is the space of functions u in H Wq’‘^{Q) satisfying conditions (j7.2l) . 

We have the following result. 


Theorem 7.3 With the notations above, we have the following trichotomy: 


i) If a > 3/2, then H, ^ 


nj- 


a) If 0 < a < 3/2 and b is non-constant, then j H, ^ 


n,D- 


Hi) If a = 3/2, then where Hq is the operator + / in fl with inter¬ 

mediate boundary conditions on dIl\W and the following boundary conditions in 
W: u = 0, Au — 7 |^ = 0, where the factor 7 is given as 


7 = 


[ \D^V\^dy = - [ bj^iAyV + AV)dy. (7.4) 

Jyx(-oo, 0 ) Jy ClyN 


Here Ay is the Laplace operator in the y-variables and the function V is Y-periodic 
in the variables y and satisfies the following microscopic problem 


A‘^V = 0, in y X (— 00 , 0 ), 

V{y,0)=b{y), on T, 

g(y,0) = 0, ony. 


(7.5) 


Remark 7.6 For completeness, we have stated in Theorem \7.3\ the three different cases 
but up to now we can only show part i) and ii). We will provide a proof of these two 
cases and will leave the proof of case a = 3/2, which is quite involved, for Section\^ 

Proof of Theorem 17.31 i) and ii). (i) Let a G]3/2,a![. It is easily verified that 
conditions (i), (ii), (iii) in Lemma 16.21 are satisfied with for e small enough, 

(ii) We prove that condition (C) is satisfied with y(H) = Wq’^{LI) and y(H,;) = 

n WQ'‘^{Vtfj- Let Kf, = VL. Notice that condition (13.2p is trivially satisfied. 
Moreover, it is easy to see that (Cl) is also satisfied. Since y(H(;) is continuously 


23 







embedded into VFQ^’^(r2e), we have that if Vf: £ with Qn^{ve) < C with C inde¬ 

pendent of e, then ||i’e||^ C for some C independent of e (see also Burenkov [TTl 
Thm. 6, p. 160]). Using Poincare inequality in the xm direction in \ fl, we easily 
get ||t'e|li 2 (n^yQ) < P(e)ll'9a:^^e|li2(Q^\f^) for some p{e) -> 0. This implies (Cl). 

We define now the extension-by-zero operator through the boundary W and 
the restriction operator to Q.. Note that is well-defined since functions in lUg’^(ll) 
vanish on W together with their gradients. With these definitions it is straightforward 
to see that conditions (C2) and (C3) (i), (ii) are satisfied. 

We now prove that condition (C3) (iii) is satisfied. Let G C 

be such that sup^^g = sup^>g ||ue||^ 2 . 2 (Q ) < oo and let v G L^(ll) be such that 

Vein ^ fo 


Note that possibly passing to a subsequence we have that v in kU^’^(ll) and 

Ve V in lU^’^(n). Moreover, proceeding as in the proof of Lemma 16.21 one can show 
that V G ITQ^’^(n). It remains to prove that Vu = 0 on W. Since u = 0 on kU we have 
that ^ = 0 on kU for alH = 1,..., — 1. Thus it suffices to show that = 0 on 

W. For this we apply Lemma 4.3 from m- To do so, for any i = 1,..., N — 1 we 
consider the vector valued function 


vP 



dVe 

dxN 


0 ,..., 0 , 


dv^ 

dxi 


(7.7) 


where the only non-zero entries are the i-th and the Wth ones. Since Ug = 0 on cffle 

then v^{x, g^{x)) = 0 for all x G W. Differentiating this last expression with respect 

(i) 

to Xi for i = 1,..., — 1 we easily get W • = 0 on T^ where v is the normal to 

the boundary. This allows to apply Lemma 4.3 from m to conclude that for any 


1 = 1,... , A^ — 1 


dv{x, 0) db{y) 
dxN dyi 


a.e. on IT X y. 


Since b is non-constant, we deduce that = 0 on lU. Thus Vu = 0 on VU hence 
condition (C) is satished. 

The proof of statement (iii) is carried out in the next section. □ 


8 Critical case ct = 3/2 (proof of Theorem 17.31 iii) 


We divide the proof into several steps organized as subsections. There are several 
important ingredients in the proof of the critical case. The first thing is to consider 
again the diffeomorphism <l>e : Og —)• 0, particularized for this situation and study 
some of its properties. This is done in Subsection 18.11 This diffeomorphism will gen¬ 
erate its pullback transformation, that we denote again by T^, which will allow us to 
transform functions defined in D to functions defined in via composition with the 
diffeomorphism 4>(;, see (18.3p below. 

With this transformation, we will consider the weak formulation of our problem 
with a test function of the type with ip a test function in D, see (18.161) below, and 
we will be able to easily pass to the limit in all terms except in a term of the type 

/ D‘^Ve ; D'^T^tfdx, 

Jwx(-e,0) 
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that requires a deeper analysis. Notice that this term carries the information of the 
oscillations of the function and the oscillations of the domain, which are coded in the 
transformation T^. Therefore, it is not surprising that this is the the most complicated 
term to analyze. We will treat this term using the unfolding operator method from 
homogenization. The definition and main ingredients of this tool particularized to our 
case are contained in Subsection 18.21 

The weak formulation of our problem and passing to the limit in all the terms, 
including the difficult one is carried out in Subsection 18.31 proving Theorem 18.401 In 
the limit problem, there is an auxiliary function h, which needs to be characterized. 
As it is customary, this is done by considering another “oscillatory test function” and 
passing to the limit appropriately. The particular calculations to characterize v are 
contained in Subsection 18.41 which concludes the proof of Theorem 17.31 Hi). 


8.1 A special transformation from to Q 


For e > 0 small, we will use the diffeomorphism ; fie —)• defined in the proof of 
Lemma 16.21 with m = 2, a = —1, g{-) = 0, and ge{-) = —e. For the convenience of the 
reader we write it here explicitly : $e(a:,a:Ar) = {x,xn — h^{x,XN)) for all {x,xn) £ 
where 

to, if — 1 < xn < —€ 


he{x,XN) 


ge{x) 


_XN±±_\^ 
g^{x)+e) ’ 


if - e< XN < ge{x ). 


( 8 . 1 ) 


We note that !>£ is in fact well-defined for any a > 1 provided e > 0 is sufficiently 
small. For this reason, although this section is devoted to the case a = 3/2, we shall 
try to keep track of a in all formulas and statements where the specific value a = 3/2 
is not required. 

The proof of the following lemma follows by straightforward computations. 


Lemma 8.2 The map is a diffeomorphism of class and there exists a constant 
c > 0 independent of e such that 


\he\ < ce", 


dh. 


dxi 


< ce 


0.-1 




dxidxj 


< cd 


0-2 


for all e > 0 sufficiently small. 


In the sequel we shall use the pull-back operator associated with Namely, 
is the operator 


T, : L2(Li) ^ T2(a) 

U !-)• tt o 


(8.3) 


Note that is a linear homeomorphism and its restrictions to the spaces LFg’^(II) and 
IT^’^(II) define linear homeomorphisms onto and respectively. In 

particular, is an isomorphism between the spaces I^(II) = IT^’^(n) H Wq’‘^{VL) and 
= IT^’^(IIe) n IlQ’"^(IIe). We also note that for any a > 1 the operator norm 
ll^e||/;(W^’^(f2) VF^’^(r2 )) uniformly bounded with respect to e, while the operator norm 
is uniformly bounded with respect to e only if a > 2. 
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Remark 8.4 As a difference from Theorem \ 7. S\ case i) and ii) in the critical case a = 
3/2 we will not he able to show condition (C). We explain here where this condition fails. 
Notice first, that is the natural candidate to show condition (C2), see Definition \3.1\ 
The other natural operator E,, : V{Qf) —^ should be the restriction operator. 

With respect to we have two different and “natural" options: K^ = Q. or K^ = Wx 
(—1,—e). For both options, condition (1321) and (Cl) hold in an easy way. Moreover, 
condition (C3) and (C2) Hi) also hold. 

The main difficulty is with condition (C2) i) or (C2) ii) depending on the choice of 
Kf:. In case = H, then (C2) i) does not hold and in case = W x (—1, —e) then 
(C2) ii) does not hold. ^45 a matter of fact it will be seen later that ||Te(/9||54/2,2(54/x(-e,o) 
does not go to 0 for most ip G 4^(0). 

Hence, it is not possible to show condition (C) for this case. 


8.2 Unfolding operator 

We will see that the limiting problem will contain an extra boundary term, which 
represents the interplay between the boundary oscillations and the boundary conditions. 
In order to identify the limiting problem in the case a = 3/2 and prove Theorem 17.31 
(hi), we shall use the unfolding operator method. In this section we recall the definition 
of the unfolding operator and some of its properties. We follow the approach of Casado- 
Diaz et al. [MlllT] and we consider an unfolding operator which is an anisotropic version 
of the classical unfolding operator discussed in Cioranescu, Damlamian and Griso m- 
We note that the well-known properties of the standard unfolding operator have to be 
slightly modified. For example, in the exact integration formula stated in Lemma 18.71 
below, an extra factor e appears in the right-hand side. Moreover, the limiting function 
V in Lemma fS. 91 below turns out to be U-periodic while, keeping in mind in the classical 
unfolding method, one would expect that the limiting function would be the sum of 
a periodic function and a polynomial of the second degree in the variables y (cf., [El 
Thm. 3.6]). 

For any k G 'L^~^ and e > 0 we consider the e-cell = ek + eY, where as above, 
the basic cell Y is given by Y = 5 )^”^- Let Iw,<i = {k G C W}. We 

set 

W,= [j Cl 

kGlw,e 

Then we recall the following 


Definition 8.5 Let u be a real-valued function defined on fl. For any e > 0 sufficiently 
small the unfolding u of u is the real-valued function defined on W^ xYx (—l/e,0) by 


u{x,y,yN) = u 



+ ey, eyw 


( 8 . 6 ) 


for all {x,y,y]y) G lU x T x (—1/e, 0), where [|] denotes the integer part of the vector 
j with respect to the unit cell Y, that is, [|] = k G if and only if x G C^ . 


We also recall the following lemma the proof of which can be carried out exactly as 
in the standard case discussed in m- 
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Lemma 8.7 (Exact integration formnla) Lat a G [—1,0) be fixed. Then 


u{x)dx = e I u{x,y)dxdy 

a,0) J WeXY x{a/e,0) 


JWeX{. 

for all u G L^(r2) and e > 0 sufficiently small. 

We denote by locO^ ^ (“Oo, 0)) the space of functions in x (—oo, 0)) 

which are also E-periodic in the first (N — 1) variables y. As customary, we shall often 
identify functions in locO^ ^ (“OO,0)) with their restrictions to E x (—oo,0). In 

the sequel we are going to consider functions in locO^ ^ (“Oo,0)) whose second 

order weak derivatives are square summable in E x (—oo,0). For this reason we hnd it 
convenient to set 


w 


2,2 

Pery 


{Y X (-00,0)) = {n G X (-00,0)) : 

ll^"^^llL 2 (yx(-oo, 0 )) < oo, Vjaj = 2} . (8.8) 


Lemma 8.9 The following statements hold: 

(i) Let Ve G kF^’^(fl) with Ijt^gjl 1 ^ 2 , 2 < M for all e > 0. Let 14 be defined by 

Veix, y) = Ve{x, y) - ^)dy - ^yVe{x, y, 0)di/ • y (8.10) 

for {x, y) G VI 4 X y X (— 1/e, 0). Then there exists v G L^(W, (E x (— 00 ,0))) 

such that 

(a) V and VyD in L‘^{W xYx (d, 0)) for any d < 0. 

(b) DyV in L^{W xY x (- 00 ,0)), for any jaj = 2, 

where it is understood that functions I4,VyI4 and DyV^ are extended by zero in 

the whole of W x Y x (—oo,0) outside their natural domain of definition x 
y X (-l/e,0). 

(a) Let V G ly^’^(n). Then 

{f^u"-^v{x,0), in L^W xY X {-1,0)) 

Proof. We start proving statement (i). It is obvious that dl"I4 = DyV^: for any 
|q:| =2. By Lemma 18.71 and the chain rule it follows that 


IWeXYx{-l/e,0) 




f3/2 


dxdy = / e|Zl“Ugpdxdy 

JWcXYx{-l/t,0) 


f |L>“ug| 2 dx < [ iD^vJ'^dx < M^, ( 8 . 11 ) 

Jwex (-1,0) Jn 


for all e > 0, hence 


is uniformly bounded with respect to e. 


L2(Wxyx(-oo,0)) 

Note that the operator dehned by fy v{y, 0)dy + Jy Vv{y, 0)dy ■ y for functions v in 
a Sobolev space of the type W‘^''^{Y x (d, 0)) with d < 0, is a projector on the space 


27 













of polynomials of the first degree in y. Thus, we can apply the Poincare-Wirtinger 
inequality and conclude that for any d < 0 there exists > 0 such that 


K 

g3/2 


L2{WxYx{d,0)) 


YvYl 

£3/2 


— ^ 

L^WxYx{d,0)) . 


£3/2 


< CdM, 

L'^(WxYx{d,0)) 


for all e > 0. A standard argument implies the existence of a real-valued function v 
defined on W xY x (—oo,0) which admits weak derivatives up to the second order 
locally in the variable y, such that v,VyV € L‘^{W xYx (d, 0)) for any d < 0, DyV G 
L‘^{W xY X (—00,0)), and such that statements (a) and (b) hold. 

It remains to prove that v is T-periodic in the variables y. Note that 


j^Vyv{x,y,0)dy = 0, 


for almost all x G W, hence it suffices to prove that is T-periodic in the variables 
y. We note that 

VyV/(x,y) _ Vv^(x,y) - fyVv^(x,y,0)dy 


p3/2 




Thus, in order conclude it is simply enough to apply the same argument in Step 3 of 
the proof of Lemma 4.3 in Casado-Diaz et al. m to the function Vn^. 

We now prove statement (ii). If n G (7°° (II), we easily can see that 


IWeXYx{-l,0) 
rO 


{Tev)\n - n(x,0) 


dxdy = 


/yx(-l,0) 


T^v — v{x, 0) 


dxdy 


E 

k€.Iw,e 


[XI 


(T,u)(e ey, eyN) - v{x, 0) 


dydxdyN 


y] / / \{Ttv){z,eyN) - v{x,[ 


kGlw, 


1 ^ -^ZidxdyN 


[ [ [ \viz,eyN - heiz,€yN)) - v{x,0)\^-^—dxdyN 

d-i Jc^ Jc^ 




< C 


E 


z - x\‘^ + \eyN - he{z,eyN)\^ -^^dxdyN 


dz 


1 - I r'k I r^k 

keiw,e 


< c [ y [ e^dxdyN < ce^, 

Jc^ 


( 8 . 12 ) 


kGlw, 


hence statement (ii) is proved for smooth functions. In the case of an arbitrary function 
V G VT^’^(n), we use an approximation argument. Namely, we consider a sequence 
Vn G C°°{Q) converging to v in IT^’^(II) as n ^ oo and we note that 

\\iTeV)\n - v{x,0)\\ < ||(Te^^)|f2 - (re^^n)loll + \\{T^Vn)\n - Vnix,0)\\ + lbn(^, 0) - n(x, 0) || 

(8.13) 

where all norms are taken in L‘^{W xYx (—1,0)). Since statement (ii) holds for 
smooth functions, the second term in the right hand-side of (I8.13P goes to zero as 
e —>• 0. Moreover, by the continuity of the trace operator, also the third term in the 
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right hand-side of (I8.13P goes to zero as e —^ 0. We now consider the first term in the 
right hand-side of (I8.13P . By Lemma 18.71 and changing variables in integrals, we have 


/ 


TfV — TfVr 


dxdy 


WeXYx{-l,0) 

= e~^ I \v o — Vn o ^t\‘^dx 

JWeX{-e,0) 

< ce 




■I-,(VK,x(-e,0)) 


V - Vnl'^dx < c||t> - Vn\\wl,2^Q) 


(8.14) 


where the last inequality is deduced by the fact that the diameter in the XAr-direction 
of the set x (—e, 0 )) is 0 (e) as e —)• 0 and that the function \v — Vn\ is bounded 

in almost all vertical lines. Since the right hand-side of (I8.14jl goes to zero as re —oo, 
we easily conclude. □ 


8.3 Weak macroscopic limiting problem 

Let /e G L?{Q.e) and / G L^(n) be such that f in L^(M'^) with the understanding 

that such functions are extended by zero outside fig and 14 respectively. Let Ue G 
17(nc) = VL^’^(ne) n WQ^’^(14e) be such that 

= fe (8.15) 


for all e > 0 small enough. By (18.1511 it follows that 11Ue|| 14 / 2 , 2 < M for all e > 0 
sufficiently small hence, possibly passing to a subsequence, there exists v G B7^’^(n) n 
VFQ^’^(n) such that Ue v in W^’^(n) and v in L?{Vl). 

Let G V{Q) = B7^’^(14) n Wg^’^(14) be a fixed test function. Since G 17(14^), 
by (j8.15p it follows that 


D^Ve : D^T^ipdx -|- / v^T^ipdx = / f^T^cpdx. 

Jue JUe 


(8.16) 


It is easy to see that 


/ v^T^ipdx / vipdx, and / f^T^ipdx / fipdx. 

Juf Jn Jn^ Jn 


(8.17) 


We now consider the first term in the left hand-side of (I8.16p . It is convenient to 
set = W X (—1, —e) so that 


/ 

Jq, 


I 

Jq 




I 

Jq\k^ 


D^Ve : D^T^ipdx = / : D^T^‘pdx+ / '■ D^T^ipdx+ / : D^T^ptdx 


-I 

Jk, 


Since D‘^T^p> = D‘^(p in we have that 


(8.18) 


/ 

JXe 


/ 

Jk, 


D^v, : D^T.ipdx = I D'^v, : D^ipdx j D^v : D'^tpdx. (8.19) 


Moreover, one can prove that 


' f2e\f2 


D^Ve : D^T.ipdx 0 . 


( 8 . 20 ) 
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Indeed, by changing variables in integrals, using the chain rule and Lemma 18.21 we 
get the following inequalities (here and in the sequel, to shorten notation we drop the 
summation symbols): 


ln^\n 

< c 

< c 

< c 


D^Ve '■ D‘^Tf:ifdx I < c 




d'^T.ip 


dxidxj 


'ne\u 










dxkdxi 


dxi dxi 


+ 


dx 




dx + ce 




dxk 

dip 




dxkdxi 


dx + ce 


-1 




dxk 

dip 


dxk 


($,(x)) 


dx. 


dxidxj 


dx 


dx 


( 8 . 21 ) 


Observe that 


^e{^e \ ^) = {{x,X]S[ — he{x,X]S[)) : X £ W, 0 < Xn < deix)} 

C {(x, zn) ■■ X £ W, -ge{x) < ZAT < 0} 

C {(x, zn) : X G ly, —e^/^6o < zn < 0}, 


( 8 . 22 ) 


where 6o = \\K')\\l°°(w)- Hence, |$e(ne \ 12)1 < ce^/2 q and therefore 


L 




d'^ip 


dxkdxi 


dx 0. 


Moreover, notice that ■§^ix,-) £ iy^’^(—1,0), a.e. x G IH and therefore 


dip 


dxi 


ix,-) 


< C 




dip 


dxi 


{xr) 


a.e. X G VL. 




Hence, the last term from (I8.2ip is analyzed as follows 


ce 


'7 




dip 


dxk 


dx < ce ^e^/^6n 


Iw 


dip 


dxk 


7>-) 


dx 


L°°(- 1 , 0 ) 


< Ce*/" 


/ 

Jw 


dip 


dxi 


{xr) 


lVi. 2 (-i,o) 


(8.23) 


We now consider the second term in the right hand-side of (j8.18p . It is convenient 
to set Qf. = We X (—e, 0) so that 

/ D^Vf_ : D^T^ipdx = / : D^T^ipdx + / : D'^T^ipdx. (8.24) 

Jn\K^ JQe Jn\{K.,uQe) 


One can prove that 


/ 


e^O 


n\{K,uQ,) 


D'^Ve : D^^T^ipdx —> 0 . 


(8.25) 
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Indeed, proceeding exactly as in ()8.21l) we get 
( [ D^Ve ; D'^T^ifdx ] 


>n\{K,uQ,) 
< c 


^,{n\{K,uQ,) 


< c 




i!>e{n\{K,UQ,) 


dxkdxi 


dx + ce 


-/ 




dip 


dxk 


dx 


dxkdxi 




(8.26) 


where now we have used the fact that the diameter of the set \ {K^ U Qe) the 
direction xn is 0(e) as e —>■ 0. Moreover, since |n\ (H^^ x (~1) 0))| —)• 0 and (/9 is a fixed 
function, we have ll<^ll^^. 2 , 2 (o\(v^ 4 x(-l,o))) ^ follows. 

It remains to analyze the first term in the right hand-side of (I8.24p . To do so we 
need the following technical lemma. 


Lemma 8.27 For all y & Y x (—1,0) and i,j = the functions h^{x,y), 

|^(x,y) and (x, y) are independent ofx. Moreover, he{x,y) = 0{e^/‘^), ^{x,y) = 

0 (eV 2 ) 

as e —)• 0 and 


1 {Ky){yN + 1 )^) 

^ dxidxj'^ ^ dyidyj 

for all i,j = 1,..., N, uniformly in y gY x (—1,0). 


(8.28) 


Proof. The independence of the functions in the statement from x is easily deduced 
by the periodicity of the function b and the definition of in (18.ip . The rest of the 
proof follows by straightforward computations and we report only those required for 
the proof of (j8.28p for the convenience of the reader. Note that 


d‘^h, 

dxidxj 


Moreover, 


.a-2 


d% 


XN + e 


+e' 


dyidyj \ge{x) +e 
db d f XN + e 


+ d 

3 


a—l 


db d 


XN + e 


dyi dxj \gt{x) + e 


dyj dxi \ge{x) + e 


+ e“fc(y)- 


Q2 


XN + e 


dxidxj \(7e(x)-|-e 


(8.29) 


d ( XN + e '\ ^ _ 36iNe Hl/Af + 1)^ _ ^(j/Af + l)^^ 

dxi \geix) + ej {e°‘-^b{y) + 1)^ {e°^-^b{y) -h l)'^ 


and 


d'^ f XN + e \ ^ _ GSiNSjNe ^(t/at -t- 1) 
dxidxj \ge{x) + e) {e°‘-^b{y) + 1)^ 

9e--HyN + imN§- + 5,wt) 

(e«-16(5) +1)! («“-'6(5) + 1)= 

{e°‘-^b{y) + 1)4 


(8.31) 
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By combining (I8.29p - (|8.3ip we easily get (j8.28p . 

We are now ready to prove the following 
Lemma 8.32 Let v € L‘^{W,Wpery{Y x (—oo,0))) be as in Lemma\8M 



D^Ve : D'^Teipdx 


€->■0 


W JYx(-1,0) 


D^vix, y) : L>^(6(y)(l + yNf)dy 


Then 

dip 

dxN 


□ 


{x,0)dx, 
(8.33) 


where Dy denotes the Hessian matrix in the variable y. 
Proof. By Lemma 18.71 and the chain rule we get 


' Qe 


; D^T^ipdx 


L 

Jw 


= e I D‘^Ve : D‘^T,ipdxdy = e“^ / 

IWeXYx(-lfi) Jw^xYx{-l,0) 

= e I ^ J ^ ^ i<^eiy))— - r^dxdy 


DyVe : DyT^ipdxdy 


+e 


-3 


j. 

/- 


M4xyx(-1,0) dyidyj dxkdxi 


dyi dyj 


IWexYx{—1,0) dy^dyj dxf^ 
We note that 


d'^Ve dip - 

a., dxdy. 


dyidyj 


d^ 


(fe) 






dyi 


- ef^, if A: = iV. 


By Lemma[8J](i), we have that II Il2,2(w? xyx(-i,o)) = 0{e^/'^) as e 0 and by 
P8.36p we have that = 0(e) as e ^ 0. Using also Lemma ISTFl once more, we get 


(8.34) 

(8.35) 

(8.36) 


.-3 


/ 


d'^Vf, d'^ip 
W4xyx(-i,o) dyidyj dxkdxi 
d^m 


($,(;/)) / dxdy 

dyi dyj 


< ce 

< ceV2 

< c 


-3/2. 


dyidyj 


L2(W^xYx{-1,0)) 




d^ip 


dxkdxi 
d^ip 


i^y)) 


< c 


dxkdxi 


L2(WExyx(-l,0)) 


dxkdxi 


L^(WeXYx{-l,0)) 

($,(x)) 


L^iQe) 


dxkdxi 


LHMQd) 


(8.37) 


Hence, the integral in p8.34p vanishes as e —>■ 0. 
We note that 


dyidyj 


0 , 


—e 


2 (92 fe, 
dxidxj ’ 


if A / iV, 
if A; = iV. 


(8.38) 


Therefore, in the sum of integrals in p8.35l) the terms corresponding to the indexes 
k ^ N vanish. It remains to analyze the term in (I8.35p with k = N. 
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By (I8.38P we rewrite such term in the form 


,-3 


lwexYx{-i,o) dyidyj dx^ 


a., dxdy 


dyidyj 


L 

JV. 


vu£xyx(-i,o) 


,-3/2 


(9^i 


dyidy. 


d(p 


.1/2, 


N 


d'^h^ 

dxjdxn 


dxdy. (8.39) 


By applying Lemma [8i9] (i) to the sequence Ug, Lemma [8/91 fii) with v replaced by 
and using (I8.‘28p we easily deduce that the integral in the right hand-side of (|8.39l) 
converges to the integral in the right-hand side of (j8.33l) as e ^ 0. □ 


Thus we have proved the following 

Theorem 8.40 Let G and f G L?‘{Q) he such that f in L‘^{Q). 

Let Ve G n WQ’'^{Lle) be the solutions to 118.15\) . Then possibly passing to a 

subsequence, there exists v G H Wq’‘^{LI) and v G L‘^{W,w‘^p‘^^^{Y x (—oo,0))) 

such that Ue v in v in L‘^{Ll) and such that statements (a) and (b) 

in Lemma \8.9\ hold, and such that 


f D‘^v : D'^if+ uipdx - [ [ Dyv{x,y) : Dy{b{y){l+ yN)^)dy-^^{x,0)dx 

Jn Jw Jyx{-i,o) 


dxM 

= f f^dx, (8.41) 
Jn 


for all if G 1L2’2(L!) n ILo’^(Ll). 


8.4 Characterization of v via a weak microscopic problem 

In this section we plan to characterize the function v dehned in Theorem 18.401 

Let ip G C°°{W xyx] —oo, 0]) be such that supp V' C CxTx [d, 0] for some compact 
set C C VL and d g] — oo, 0[ and such that f^{x, y, 0) = 0 for all (x, y) £ W xY. Assume 
also that if is T-periodic in the second {N — l)-variables y. We set 

ifeix) = e^ip (x,^,^^ , (8.42) 

for all e > 0, X G VLx] — oo, 0]. Note that for e sufficiently small we have that supp V'e C 
n and 'i/’e G V{Q), hence belongs to V{Qe) and can be used as test function in the 
weak formulation of our problem in fig, that is. 


/ : D'^T^fj^dx + / ufT^fj^dx = / ffT^^if^^dx. 

'n^ Jn^ 


I 

Jn. 


By the presence of the factor e 2 in (I8.42h , it is easy to see that 


/ 

Jn, 


e^O 


vfTf^'ifedx —> 0, and / feT^^ipf^dx 


e^-O 


0 . 


(8.43) 


(8.44) 


We now consider the hrst term in the left hand-side of (|8.43l) and we write it in the 
form 


: H^Te^edx = / D^Ve '■ D'^T^^if^dx + / D^Ve : Il^Te'i/jedx (8.45) 

J SI 
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With very similar arguments as the ones used to show ()8.20p we can also prove that 


We can show now the following 


[ D\ : D‘^T,iP,dx ^ 0 . 

Jn, 


(8.46) 


Lemma 8.47 is as in (8.4-2) andve, v are the functions from Theorem\8.40[ then 


/ 


D^Ve : D^T^fj^dx 


Dlv{x,y) : Dy'ilj{x,y)dxdy. 


JWxY x(—oo,0) 

Proof. First of all we note that by the periodicity of if we have that 


T^ifeix^y) = 

We also note explicitly that 


X 


L e J 


+ f^y,y,yN - e ^K[e 


X 

L e J 


+ ey, (-yN 


(8.48) 


(8.49) 


h, (e 


X 
e J 


+ ey, eyAf = 


(2''^b(y){yN+lf If _ 1 < / n 

(el/26(5) + l)3 > II J- < yAT < U, 

0, if - 1/e <yN < -1, 


hence 




^e ^e — + ey, eyAT^^ 0, uniformly on W x y x] — oo, 0], (8.50) 


for all |a| < 2 . 

Since -0 is smooth and has compact support, it is Lipschitz continuous together with 
its derivatives and it easily follows that 


(e 


X 

L e J 


+ (^y,y,yN - e ^hAe 


X 
e J 


+ ey,eyN)) - Df^'tf{x,y) 

)) L 2 (^^xyx]-oo, 0 [) 

(8.51) 

for any |/3| < 2. In fact, the square of the norm in (|8.51l) can be estimated by 


dxdy 


e^O 


f 


'X 


2 

e ^he (e 

'X 

\ 

L 

e 

— 

ey — X 

+ 

— 

+ ey, eyN 

JWeXYx]-d,0[) 





\ 

. e. 

/ 


which is clearly 0(e) as e ^ 0 . 

By combining (I8.49|) - (|8.5ip and using the chain rule, we get that 

e-y^D()T,'ilj,(x,y) ^ D(f'if(x,y) 

in L‘^{W xYx]— oo, 0 [), for all I 7 I = 2. 

By Lemmas 18.7118.91 and (I8.52P , we conclude that 


(8.52) 


/ 


: D‘^Tf:iff:dx 


ty£x]-i,o[ 


>WeXYx]-l/e,0[ 


£3/2 


D'^.T.v, _ 


£3/2 


-dxdy 


'W xY x(—oo,0) 


Dyv(x,y) : Dy'il;{x,y)dxdy 


□ 
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Theorem 8.53 (Characterization of v via a two scale weak problem) Let v G 

X (—00,0))) he the function from Theorem \8.40\ Then 


/ Dlv{x, y) : y)dxdy = 0 , (8.54) 

JlVxVx(-oo,0) 

for all if G L'^(W,w'^p1j.^(Y x (—oo,0))) such that 'if{x,y,0) = 0 on W xY. Moreover, 
for any i = 1,..., N — 1, we have 

y, 0) = ^{y)-^{x, 0), on IT x T . (8.55) 

dyi dyi dxN 

Proof. The proof of (18.541) for smooth test functions if follows by passing to the 
limit in equation (I8.43P and combining (|8.43l) - (|8.48p . The general case involving test 
functions if G L^{W, {Y x (—oo, 0))) follows by an approximation argument which 

we skip for brevity (we only mention that in order to preserve the boundary condition 
at yN = 0, one can first extend a given test function if by setting if{x,y,—yN) = 
—if{x,y,yiy) and then using convolution). 

We now prove (I8.55p . Let = ^0,..., 0, ..., 0, be the vector as in 

_ 

(j7.7j) for any i = 1 ,..., — 1. As already noted, K • = 0 on Lg where u is the 

normal to the boundary. We note that by Lemma 18.91 we have that 

d'^v 

Ve dyj I dxi j dyjdyi 


in L‘^{W X Tx] — oo,0[). Thus, it suffices to apply Casado Diaz et al [TFl Lemma 4.3] 
to conclude. □ 

We plan now to describe function t) in a more explicit way by separating the variables 
X and y and providing a classical formulation of the microscopic problem (|8.54p . To do 
so, we shall also need to perform some calculations based on standard integration by 
parts on domains of the type Y x (d, 0) with d < 0. In essence, the computations are the 
same computations required to prove a known Green-type formula for the biharmonic 
operator. For the sake of clarity, we state such formula and we provide a short proof in 
the classical setting. We recall that, given a smooth open set U and a smooth vector held 
F : dU the tangential divergence of F is dehned by divg[/T = divT — (VT -u)!/ 

where it is meant that the vector held F is extended smoothly in a neighborhood of 
dU. In the following statement, we shall also denote by Fgu the tangential component 
of a vetor held F as above, dehned by Fgu = F — {F ■ v)iy. 

Lemma 8.56 (Biharmonic Green Formnla) Let U be a bounded open set in 
of class . Let f G C^{U) and ip G C^{U). Then 

I Dy:D\<ly=Uy^ay-l 

Ju Ju JdU dv Jgu ov OV 

+ f • ^)du ■ Yguifda. (8.57) 

JdU 
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Moreover, if Vt is also of class then 


i OV : D^fdy = (div8„((DV) ' -)» + Vda. 

(8.58) 


Proof. Standard integration by parts gives 



D^f : D^^dy 



d^f d'^ip 
dyidyjdyidyj ^ 



d^f dip , f d^f dp , 

O 2 a / Q a 

dyfdyj dyj Jqu dyidyj dyj 


lu 


AfApdy - 


Idu dv 


'du 


{D^f • n)Vpda, 


(8.59) 


where summation symbols have been dropped. Formula (18.571) simply follows by the 
standard Green formula and by decomposing \7p as Vg[/<^ + in (I8.59p . Finally, 
formula (j8.58p follows by applying the Tangential Green Formula (see, e.g., Delfour and 
Zolesio [23l §5.5]) to the last integral in the right-hand side of (I8.57p . □ 


We also need the following lemma where by Wp^^^{Y x (d, 0)) we denote the space 
of functions in x (d, 0)) T-periodic in the first {N — l)-variables y and such 

that all derivatives up to the fourth order are square summable in P x (d, 0 ). 

Lemma 8.60 There exists V € ^ (“Oo,0)) satisfying the equation 


D^V 


>Y x(—oo,0) 


D'^fjdy = 0 , 


(8.61) 


for all if G WplryO^ ^ (-oo,0)) with 'if{y,0) 


0 on Y, and the boundary condition 


V{y,0)=b{y), on y. 


(8.62) 


Function V is unique up to the sum of a monomial of the type ayn where a is any real 
number. 

Moreover, V is of class Wp'^.^^{Y x (d,0)) for any d < 0, A’^V = 0 and it satisfies 
the boundary condition 

d'^V 

-^(y,0) = 0, on y. (8.63) 

Proof. In order to prove the existence of V, we use a standard direct method in 
the Calculus of Variations. Let A = {V G w^Pery(V x (—oo,0)) : y(y,0) = h{y) on y}. 
Obviously, It is clear that a minimizer (if it exists) for the problem 

inf [ |L>Vpdy (8.64) 

Jyx(-oo,0) 

is a solution to problem (18.611) and satisfies the required boundary condition (see also 
Lemma I8.56P . Thus it suffices to prove that there exists a minimizer in (I8.64p . Let 
G A, n G N be a minimizing sequence for (j8.64p . Since this sequence is bounded in 
X (— 00 , 0 )) there exists V G w'peryO^ ^ (“Oo,0)) such that D^Vn converges 
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weakly to D^V in L‘^{Y x (—oo,0)). Since the trace operator is compact we have 
that Vn{y,0) converges strongly to V{y,0) in L'^{Y), hence V{y,0) = b{y) on Y and 
V E A. We now prove that V is in fact a minimizer. By the elementary inequality 
|ap > |/3p + 2/3 ■ {a — f3) valid for all vectors a, f3 in any Euclidean space, we get 

[ \D^Vn\‘^dy> [ \D^V\‘^dy + 2 [ : {D^Vn - D'^V)dy. 

Jyx{—oo,0) Jyx{—oo,0) Jyx{—oo,0) 

Passing to the limit in the previous inequality, using the weak convergence and the fact 
that the sequence is minimizing in (I8.64p . we find out that 

inf [ \DW\‘^dy= [ \D'^V\'^dy, 

JVx(-cx),0) Jyx{-oo,0) 


hence E is a minimizer. 

Uniqueness is easily proved by observing that if U € '^Pery^ ^ (“00,0)) satisfies 
equation (I8.61jl and the boundary condition V (y, 0) = 0 on Y, then V itself can be tested 
in (I8.6ip . and it easily follows that D^V = 0. This implies that U is a polynomial of 
the first degree of the type ^he periodicity in y we get aj = 0 for all 

j = 1,..., N — 1 and the proof of the first part of the statement is complete. 

Regularity of V is standard. The rest of the proof follows by using in the weak 
formulation (|8.61|) test functions i/j as in the statement, with bounded support in the 
yAT-direction and using formula (18.571) . In fact, using such test functions 'll), we get that 
the boundary terms corresponding to the subset dY x (—oo, 0) of the boundary cancel 
out because of periodicity, hence 


/ 


Y X (— 00 , 0 ) 


D^V : D^'il)dy = 


'Y x(—oo,0) 


{^^V/'ipdy + ^ 


d‘^V di/) 
Y du'^ du 


da = 0. 


By the previous formula and the arbitrary choice of t/: we deduce that V is biharmonic 
and satisfies the boundary condition (I8.63p . □ 


Lemma 8.65 Let V be as in Lemm,a \8.6(A Then 

[ D^V:D^{b{y){l+yNf)dy=[ \D^V\^dy = - [ bJ-{AyV+AV)dy. 

Jyx{-1,0) JWx]-oo,0[ Jy <yyN 

( 8 . 66 ) 


Proof. Let i/) be the real-valued function defined on Tx] — oo, 0[by 


^(y) 


b(y)(l + yAr)^ if (y,yAf) G u X [-1,0), 

0, if (y,yAr) G Ex]-00,-1). 


Clearly, G w’pIrY (^ ^ (“O®; 0))) i’iy^ 0) = U(y, 0). The proof of the hrst equality 
in (j8.66p then easily follows by using tp = V — ijj as a test function in (I8.61h . 

In order to prove the second equality in (I8.66p . we use formula (I8.57P with U = 
Y X (—1,0), f = V and ip{y,yi\f) = b{y){l + yN)^- First of all, we note that 
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on the part of the boundary given by Y x {0}. Then by Lemma 18.601 and formula 
(j8.57l) . by exploiting the periodicity of V and b which allows to get rid of the boundary 
terms corresponding to the subset of the boundary dY x (—1,0), we have that 

[ D^V : D\b{m + yNf)dy = [ ^y( ' V,-6(y)dy - [ b^dy. 
JYxi-lfl) Jy V (^VN J ' Jy (JUN 

Integrating by parts the right-hand side of the previous equality and using again the 
periodicity of the functions we conclude. □ 


Finally, we can prove the following 


Theorem 8.67 Let V be as in Lew,ma \8. 6 (A Let v,v be the functions defined in Theo- 
rem \8.40\ Then v{x, y) = Viy)-^^{x, 0)+a{x)yN for almost all {x,y) gWxY x]—oo, 0] 
where a(-) G L^(VF) is a function depending only on x. 

Moreover, for the strange term in hSfil] ) we have 


'w Jyx{-i,o) 


Dyv{x,y) : Dl{b{y){l + yNf)dy-^{x,0)dx 


|T>' 


/rvx]-oo,0[ Jw dxN 

OXJM OXJM 


V\‘^dy f (x,0)dx ^^ (x,0)dx 
Jw 


dxN 
dxN 


( 8 . 68 ) 


where 7 is defined in (B- 


Proof. Recall that function v satisfies problem (I8.54p and the boundary condition 
(I8.55P . By proceeding exactly as in the proof of Lemma 18.601 one can easily see that 
such function v is unique up to the sum of a function of the type a{x)yN as in the 
statement. The proof then follows simply by observing that a function of the type 
l/(y)^^(x,0) as in the statement satisfies problem (I8.54p and the boundary condition 
(I8.55p . Equality (j8.68p follows by ()8.66l) . □ 
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